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Abstract 

We study mutual dissociation of heavy nuclei in peripheral collisions at 
ultrarelativistic energies. Earlier this process was proposed for beam lumi- 
nosity monitoring via simultaneous registration of forward and backward 
neutrons in zero degree calorimeters at Relativistic Heavy Ion Collider. 
Electromagnetic dissociation of heavy ions is considered in the framework 
of the Weizsacker- Williams method and simulated by the RELDIS code. 
Photoneutron cross sections measured in different experiments and calcu- 
lated by the GNASH code are used as input for the calculations of dissocia- 
tion cross sections. The difference in results obtained with different inputs 
provides a realistic estimation for the systematic uncertainty of the lumi- 
nosity monitoring method. Contribution to simultaneous neutron emission 
due to grazing nuclear interactions is calculated within the abrasion model. 
Good description of CERN SPS experimental data on Au and Pb dissocia- 
tion gives confidence in predictive power of the model for AuAu and PbPb 
collisions at RHIC and LHC. 
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1 Introduction 



The study of a new form of strongly interacting matter, the so-called quark-gluon 
plasma, is at the core of current and future experimental programs at Relativistic 
Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) |l|] and 
Large Hadron Collider (LHC) at CERN Although colliders give well-known 
advantages compared to the fixed target experiments, the kinematics of ultra- 
relativistic heavy-ion collisions at colliders creates certain complications in the 
beam monitoring as well as in the identification of collision events. 

Due to the geometrical factor 27r6, where b is the impact parameter, the 
number of central nuclear collisions {b ^ 0) is relatively very small in the whole 
set of the collisions with nuclear overlap, b < Ri + R2 {Ri and R2 are the 
nuclear radii). Moreover, in peripheral collisions without direct overlap of nuclear 
densities, b > Ri + R2, one or both nuclei may be disintegrated by the long-range 
electromagnetic forces. This process of Electromagnetic Dissociation (ED) is a 
well-known phenomenon 0, ^. The properties of central and peripheral collisions 
are very different and then should be studied separately. The ED events are less 
violent than the collisions with the participation of strong interactions. Namely, 
the average particle multiplicities are essentially lower ^ and the main part of 
nucleons and mesons is produced in projectile and target fragmentation regions, 
very far from the mid-rapidity region. 

Calculations show ^, ^ that the ED cross section in collisions of heavy 
nuclei at RHIC and LHC by far exceeds the dissociation cross section due to the 
direct nuclear overlap. In AuAu and PbPb collisions at such energies many neu- 
trons can be produced in the ED process . Among other interesting phenomena 
one may expect a complete disintegration of nuclei induced by the electromag- 
netic fields of collision partners This phenomenon is very well known in 
nuclear reactions under the name of "multifragmentation" ^j. 

Several operational problems of heavy-ion colliders are connected with the 
high rate of the ED process. On the one hand, the ED process reduces the 
lifetimes of heavy ion beams in colliders as compared with the proton-proton 
accelerator mode ^, . On the other hand, the process of simultaneous neutron 
emission from the collision partners, where the ED process plays a dominant role. 



can be useful for the luminosity monitoring [10], pJ]] . 

The luminosity monitoring method based on mutual dissociation has several 
advantages |^ |TD|, |rT|. In particular, the beam-residual-gas interaction events 



can be strongly suppressed in favour of the beam-beam events by the condition 
that a pair of neutrons should be detected in coincidence by each arm of the 
calorimeter. The cross section of mutual neutron emission can be calculated 
in the framework of conventional theoretical models designed for describing the 
heavy ion disintegration in peripheral collisions. Corresponding nuclear data, 
especially photoneutron emission cross sections, may be used as numerical input 
for such calculations. Therefore, the neutron counting rates in zero degree (very 



3 



forward) calorimeters may provide an accurate measure for the heavy-ion colhder 
luminosity. 

In the present paper the neutron emission in peripheral collisions of ultrarel- 
ativistic heavy ions is considered with the aim of providing the theoretical basis 



for the luminosity monitoring method proposed in Refs. |iy, |TT|. The uncer- 
tainties in results originating from uncertainties in input nuclear data and in the 
theoretical model itself are carefully examined. A brief review of corresponding 
photonuclear data is given with special attention to the publications describing 
data evaluation and re-measurement. Model predictions for the Au and Pb frag- 
mentation cross sections are compared with recent experimental data obtained in 
fixed target experiments at CERN SPS with the highest energies available thus 
far. This serves as an important test before extrapolating our methods to the 
RHIC and LHC energies. 



2 Equivalent photon approach to simultaneous 
electromagnetic dissociation 

2.1 First order dissociation processes 

The electromagnetic excitation of one of the collision partners, A2, followed by 
its dissociation is schematically shown in Fig. |l]. In such a process another part- 
ner, Ai, emits a photon, but remains in the ground state without any nuclear 
excitation. Besides this "classical" process, one can consider a non-classical pro- 
cess where the emission of a photon is accompanied by the nuclear excitation 
(see Fig. |^), particularly, the giant resonance excitation. Such lowest-order con- 
tribution to the simultaneous (mutual) excitation of the nuclei Ai and A2 was 
considered in Refs. [|12|, |13|. Also the correction to photon-photon luminosity 



function due to the inelastic photon emission was considered in Ref. for 77 
fusion reactions. 

As shown in Refs. [0, |l3l (see also the discussion in Ref. @]), the lowest 
order process of simultaneous excitation of the collision partners has a small 
cross section. For the cases of interest, i.e. for AuAu and PbPb collisions, the 
cross sections for the simultaneous dipole-dipole excitation of such nuclei are 0.49 
and 0.54 mb, respectively |]12|. Using a rough estimation of Ref. |Q, 10~^A'^ mb, 



one can get for the same nuclei 0.39 and 0.43 mb, respectively. From the following 
discussion one will see that for heavy nuclei these first order contributions (Fig. |^) 
are negligible compared to the second order ones (Fig. The latter leading order 
process with exchange of two photons is a classically allowed mechanism. It is 
considered in the next section where the formalism previously used in Ref. is 
extended to the case of mutual excitation. 
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2.2 Second order dissociation processes 

Let us consider a collision of heavy ultrarelativistic nuclei at impact parameter 
b > R1+R2. The masses and charges of these nuclei are denoted as Ai, Z\ and A^, 
Z2, respectively. Hereafter the case of equal nuclei {A\ = A2 = A, Z\ = Z2 = Z 
and Ri = R2 = R) is investigated. Nevertheless, in some cases the indexes are 
used to show explicitly which of the collision partners emits or absorbs photons. 



According to the Weizsacker- Williams (WW) method JTSl, the impact of the 
Lorentz-boosted Coulomb field of the nucleus Ai on A2 is treated as the absorp- 
tion of an equivalent photon by the nucleus A2 (see Fig. In the rest frame 
of this nucleus the spectrum of virtual photons from the collision partner Ai at 
impact parameter b is expressed as: 

N.AEub) = ^a4^(A-?W + :^A-(x)). (1) 

Here a is the fine structure constant, x = Eib/{'yf3hc) is an argument of the 
modified Bessel functions of zero and first orders, Kq and Ki, j3 = v/c and 
7 = (1 — is the Lorentz factor of the moving charge Zi. If the Lorentz 

factor of each heavy-ion beam is 'jbeam, then 7 = 27^^^^ — 1 for the case of collider. 
Hereafter the natural units are used with h = c = 1. 

The mean number of photons absorbed by the nucleus A2 in the collision at 
impact parameter b is defined by: 

mAM= J Nz,{E^,byA,{E^)dE^. (2) 

where the appropriate total photoabsorption cross section, aA2{Ei) is used. For 
Emin one usually takes the neutron emission threshold, while the upper limit 
of integration is E^ax ~ 7/-R. We assume that the probability of multiphoton 
absorption is given by the Poisson distribution with the mean multiplicity (^) 
defined by Eq. (1). 

Following Refs. R ITB, O , we express the cross section for the electromagnetic 



dissociation of one of the nuclei due to the absorption of a single photon (Fig. |l|) 
leading to a certain dissociation channel i as: 

00 

af^(z) = 27r J bdbPA,{b), (3) 

bo 

where the probability of dissociation at impact parameter b is defined as: 

F, 

PaM = e-"^''^'^ J dE,Nz,{E,,b)aA,mfA,{E,,i), (4) 
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and fA2iEi,i) is the branching ratio for the considered channel i in the absorption 
of a photon with energy Ei on nucleus A2. The choice of a critical impact param- 
eter be, which separates the domains of nuclear and electromagnetic interactions, 



will be discussed in Sec. 4.4. 



Let us turn now to the mutual dissociation process shown in Fig. |^. The cor- 
responding graph may be constructed from two graphs of the single dissociation 
by interchanging the roles of "emitter" and "absorber" at the secondary photon 
exchange. Several assumptions have to be made to obtain an expression for the 
mutual dissociation cross section. 

Firstly, we suppose that the emitted photon with energy Ei < Emax does not 
change essentially the total energy, Ea = 'J Ma, of the emitting nucleus, where 
Ma is the nuclear mass. This can be justified by estimating the ratio 

which is close to 10~^ for heavy nuclei. Therefore, the kinematical conditions 
for the secondary photon exchange are very similar to those for the primary one, 
and there are no correlations between the energies of the primary and secondary 
photons. El and £'2. In other words, the primary and secondary photon exchanges 
may be considered as independent processes even if they take place in the same 
collision during a short-term overlap of the Lorentz-contracted Coulomb fields of 
the colliding nuclei. Secondly, the equivalent photon spectrum from the excited 
nucleus, A2 in the notations of Fig. 0, is the same as the spectrum from the nucleus 
in its ground state, A2. This follows from the fact that at ultrarelativistic energies 
the collision time is much shorter then the characteristic deexcitation time during 
which a nucleus changes its initial charge via proton emission or fission. 

Following these assumptions, one can express the cross section for the mutual 
dissociation of nuclei Ai and A2 (Fig. ^ to channels i and j, respectively, as: 

00 

I j) = 2nJ bdbPA,{b)PA,{b). (6) 

fee 

Substituting Eq. (^) for each of the nuclei and changing the order of integration 
one obtains: 

^max Emax 



(^nF{^\j)= j j dEidE2ArUEi,E2)aA,{Ei)crAAE2)fA,{Ei,t)fAAE2,j), 

(7) 

where the spectral function Mm{Ei, E2) for mutual dissociation is introduced: 

Mm{Ei, E2) = 2ttJ bdbe-'^'-'^Nz.iEi, b)Nz,{E2, b). (8) 
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Conditions A\ = = A and Z\ = Zi = Z were used in Eqs. (0) and (Q), 
as it is usually in heavy-ion colliders, and therefore rriA^i^ = mA2{b) = m{h). 
Nevertheless, the dissociation channels i and j may be different for each of the 
nuclei even in such case. 

Several remarks may be made concerning Eqs. (0) and (^. Compared to 
another process of the second order dissociation of a single nucleus (see Fig. ^) 
some points of similarity may be found. Indeed, expressions given in Ref. 
for the corresponding cross section o"|'^(z) of the second order process are as 
following: 

(^2''i^= I I dEidE2U2{E^,E2)aA,{E^)crAAE2)fA,{E^.E2,i), (9) 

OO 

Ar2{E,,E2) = nj bdbe-"''^'^Nz,{E,,b)Nz,{E2,b). (10) 

be 

However, there is an important difference in the definitions of branching ratios 
fA2{Ei,i), fAi{E2,j) compared with fA2{Ei,E2,i), since the former are for the 
absorption of two photons by two different nuclei leading to certain dissociation 
channels i and j, while the latter is for the absorption of two photons by a single 
nucleus leading to a channel i. Another difference is due to an additional factor 
of 2e-'^('') in Eq. (^ compared with Eq. (1). It comes from the fact that Eq. (^) 
contains the product of the Poisson probabihties for the single photon absorption 
for the collision with impact parameter b: 

PaMPaM = m2(6)e-2™('') = 2e-™(^)pi'^(6), (11) 

while PA\b) = m^(6)e^™*^^''/2 is the Poisson probability for the double photon 
absorption |16|. Results for partial cross sections of mutual dissociation will 
be presented and discussed in Sec. ^. 



2.3 Next-to-leading-order processes of mutual dissocia- 
tion 



The second order process, which is shown in Fig. |] and considered in Sec. ^ 



is the leading order mechanism of mutual electromagnetic dissociation. This is 
confirmed by the calculations of the corresponding total cross sections for the 
leading order (LO), Fig. ^ and next-to-leading-order (NLO) processes. Figs. ^ 
and p. 

Following the assumption that the probability of multiphoton absorption is 
given by the Poisson distribution with the mean multiplicity mA{b) of Eq.(0), one 
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has for the LO process shown in Fig. |^: 



CT^^(LO) = 271 J hdhm\{h)e-^'^^^^\ (12) 

6c 

Here the case of equal masses and charges of colhsion partners is considered, i.e. 
Ai = A2 = A and Zi = Z2 = Z. 

For the NLO process with exchange of three photons (NLO12), Fig. ^, the 
total cross section is given by: 

a;^^(NL0i2) = 27, j hdh'^^e-^^^^'\ (13) 

be 

A complementary process (NLO21) with the excitation of another nucleus, A2, 
via double photon absorption is equally possible and has the same cross section. 

Another NLO process of mutual dissociation is due to exchange of four pho- 
tons (NLO22), Fig. and its cross section can be written as: 

a^^(NL022) = 27iJbdb^^^e-^^^^'\ (14) 

be 

Finally, one can calculate the sum of all contributions using the prescription 
of Ref. II: 

a^^(tot) = 271 J bdb(^l - e-'"-^^^)) . (15) 

6c 

Since for each of the ions the probability of collision without photon exchange is 
equal to e~'^'^^^\ Eq.(p^ is evident. 

Calculations of cr^-^(tot), cr^^(LO), cr^^(NLO) and cross sections for specific 
dissociation channels, cr^ii \ j), with and without NLO corrections were per- 
formed by the modified RELDIS code, which contains now a special simulation 
mode for the mutual electromagnetic dissociation process. Our results for the 
total dissociation cross sections in electromagnetic and nuclear interactions are 
given in Tab. |l[ These cross sections values are in good agreement with results 
of other authors. 

The total cross sections for mutual electromagnetic dissociation given in Tab. |l] 
are much lower than the cross sections for single dissociation. However, even the 
former values are found to be comparable to the total nuclear dissociation cross 
sections, see Tab. |I|. 

As one can see, the ratios between the first and second order processes are 
very different for single and mutual dissociation. The first order dissociation 
process (with exchange of a single photon) can be safely neglected in considering 
the mutual dissociation of heavy nuclei at ultrarelativistic energies. 
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Tab. |T] contains also the cross sections for LO process, a^^{LO). As one can 
see, the LO mechanism gives ~ 70% and ~ 60% of the (j^^(tot) at RHIC and 
LHC energies, respectively. The sum of the NLO contributions to the total cross 
section gives additional ~ 25-27%. Therefore, at RHIC energies, for example, the 
remaining contribution of ~ 5% (or ~ 0.2 b) is due to exotic triple excitations. 

Finally, and Afm written for the LO mutual dissociation process, Eqs. (|^) 
and (1), can be easily generalized to the case of different NLO processes. However, 
the resulting expressions are lengthy and we do not give it here for brevity's sake. 

3 Abrasion model for mutual dissociation in nu- 
clear collisions 

Several nucleons can be abraded from collision partners in grazing nuclear colli- 
sions. We are interested in a situation when only few nucleons are removed. This 
is the case when nuclear densities overlap weakly and mainly nuclear periphery 
is involved in the interaction. 

The cross section for the abrasion of a nucleons from the projectile [Ai, Zi) 
in a collision with the target {A2, Z2) may be derived from the Glauber multiple 
scattering theory |T^: 



a""'=(a) = ( ^] x27ij hdh{l - P{b)yP{bf'~^. (16) 



— * 

Here P{b) is calculated as the overlap of projectile, pi{r), and target p2{r) den- 
sities in the collision with impact parameter b: 

P{b) = j (fsDi{s)exp(-A2aMND2{s + b)), (17) 

where the nuclear thickness functions, 



-Dl,2(s) = j rfzpi,2(s, Z) 



are introduced. In our calculations the nuclear density functions are approxi- 
mated by Fermi functions: 

Pi.2(r) = ^° ./3, . (19) 

l + exp(^) 

where Tq is a parameter which defines the nuclear half-density radius, Ro = 
To X A^/^, and d = 0.54 fm is the diffuseness parameter. The choice of these and 
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other important parameters of the model, the integration cutoff parameter, b, 



and the total nucleon-nucleon cross section, (t^n, is discussed in Sec. |44 . 

The above expressions determine only the number of nucleons removed from 
the projectile and do not specify how many protons or neutrons were knocked 
out. Further assumptions are needed to determine the charge-to-mass ratio of the 
residual nucleus and hence the numbers of protons, 2, and neutrons, n, abraded 



from the initial nucleus (see also the discussion in Ref. PQ])- 

In the present work we use the so-called hypergeometrical model |1T9| , pO|] , 
assuming that each removed projectile nucleon has a Ni/Ai probability (A^^i = 
Ai — Zi) of being a neutron: 

a'^"^(n,^) = ^7 \"^ ff"""(a). (20) 



In other words this means that there is no correlation at all between the proton 
and neutron distributions and the abrasion process removes protons and neutrons 
from the projectile nucleus in a random way. 

Several physical processes, which might be important in heavy-ion collisions, 
were neglected in this model. The excited residual nucleus created due to the 
abrasion process should undergo its de-excitation on the second ablation step. 
On this step more neutrons may be emitted via evaporation. However, as it was 



noticed in Ref. |T^, the excitation energies obtained on abrasion step due to 
removal of one or two nucleons are generally not sufficient for intensive particle 
evaporation. Therefore, for the cases of interest, i.e. In and 2n dissociation 
channels, the ablation step can be neglected. 

The abrasion of nucleons from projectile and target proceeds via high-energy 
collisions between nucleons. Nucleon-antinucleon pairs may be created in such 
interactions and neutrons may be presented in these pairs. However, as one can 



find in a compilation [^], even at high energies ~ 100 AGeV, the relative rate of 
such pair production is not so high, ~ 5%. Because of this, we do not consider 
such processes in calculations of the neutron emission cross sections. 

Knocked-out nucleons can also suffer a final state interaction with specta- 
tors |19[. We believe that this process is less important at high energies com- 



pared with intermediate energies of ~ 0.1 — 1 AGeV. For the latter case the 
escape probability is estimated to be Peso ~ 0.5 — 0.75 for peripheral collisions of 



heavy nuclei |22[- The momenta of recoil nucleons may be comparable with the 
Fermi momentum of intranuclear nucleons and their angular distribution is very 
wide so that they can be easily captured by one of the spectators. The situation 
is different at high energies, where the transverse momenta of collided nucleons 
are typically large, of order of 0.5 GeV, and therefore their subsequent capture is 
less probable. Other effects like a finite hadronization length may further reduce 
the secondary interaction probability. Therefore, we assume that in peripheral 
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collisions at RHIC and LHC the probability for each of the collided nucleons to 
escape the residual nuclei is close to unity, Pesc ~ 1- 

As we will show in Sec. ^ the above-mentioned simplifications do not lead to 
noticeable disagreements with experimental data when the removal of one, two 
or three nucleons is considered. It means that either above mentioned physical 
effects are negligible, or they compensate each other in peripheral nuclear colli- 
sions with removal of only few nucleons. However, the predictions of the present 
model for more central collisions with removal of many nucleons should be taken 
with care. 

This simple abrasion model can be easily extended to the case of mutual 
dissociation. The cross section for the removal of ni neutrons and zi protons 
from the projectile (A^^i, Zi) simultaneously with the removal of n2 neutrons and 
Z2 protons from the target (A^2; ^2)5 (-^2 = A2 — Z2) may be written as: 



ar (ni, z, I n2, ^2) = x a-(a) x (21) 



(5 


)G) 


( 


*) 


1 



Since the number of nucleon-nucleon collisions in such a process is assumed to be 
equal to a, the condition zi + rii = a = Z2 + n2 holds. 

Using this condition for the process of mutual dissociation with given numbers 
of neutrons, rii and n2, removed from the projectile and target nuclei, respectively, 
one has finally: 

I n2) =E^""^(^i.^i) ^"^?r:i. r • (22) 



\ zx+n-i ) 



Here the cross section for single abrasion process, (T"'"'^(ni, zi), given by Eqs. 
and fEDD was used. 



4 Input data for heavy ion dissociation calcula- 
tions 

As shown in Sec. |^, the photonuclear cross sections are used as input data in 
calculations of the electromagnetic dissociation cross sections. This is verified by 
the coherent nature of the photon emission by the collision partner as a whole. 
Since these photons represent the Lorentz-boosted Coulomb fields of original 
nuclei, their virtuality is very small, < l/R"^, i.e. these photons are almost real. 
Therefore, one can use the photonuclear reaction data obtained in experiments 
with real monoenergetic photons and apply theoretical models describing such 
photonuclear reactions. 

The accuracy of the mutual dissociation cross section calculation depends 
heavily on the quality of the data and parameters used as input. As we found. 
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for example, the mutual dissociation cross section is more sensitive to the proper 
choice of the critical impact parameter, be, and to the input photonuclear cross 
sections than the single dissociation cross section. The input data and model 
parameters are discussed in detail in the following sections. 

4.1 Photoneutron cross sections measured in experiments 
with real photons 

Over the years, the photoneutron cross sections for different nuclei have been 
measured with monoenergetic photons at Saclay [25[| and Livermore 



Data on different cross sections obtained in these and other laboratories were 



collected in compilations of Refs. |27, 28 



Concerning the nuclei of interest, ^^^Au and ^''^Pb, the detailed data were 
obtained mainly for (7, n) and (7, 2n) reactions, while less detailed data exist 
for (7, 3n) and (7,4n) reactions, see Refs. ||2^, The measurements were 



performed in the photon energy region 6 < E'^ < 35 MeV, where the excitation 
of giant resonances plays a dominant role, see Figs. |^ and |[ At such energies, 
the emission of charged particles {p, d,^ He,^ He) is suppressed by a high Coulomb 
barrier of heavy nuclei. Therefore, the sum of partial cross sections for all neutron 
multiphcities: (7(7, n) + (7(7, 2n) + cr(7, 3n) + cr(7, 4n) nearly coincides with the 
total photoabsorption cross section. Each of the inclusive cross sections (7(7, m) 
includes a small contribution from the channels with charged particles, (7,mp), 
(7,m2p),... At the same time such channels as (7,p), (7, 2p) were neglected at 
all. According to Ref. [^, this leads to a small systematic error ~ 3 — 5% in 
the total photoabsorption cross section measured at low energies via the neutron 
detection. 

Above the giant resonance region, at 35 < -E^ < 140 MeV, the quasideuteron 
mechanism of photon absorption dominates. Only average characteristics of pho- 
ton absorption by ^^'^Pb were measured in Saclay [^, 25 1 in this energy region. 
Neutron yields, J2i>iic!'{'y,i^), and the cross sections for emission of at least j 
neutrons, Y.i>ji(^{l were obtained in addition to the mean value and the 
width of the neutron multiplicity distribution. 

To the best of our knowledge there are no direct measurements of neutron 
emission cross sections and multiplicities at > 30 MeV in photoabsorption 
on ^^^Au. The only attempt to deduce the average photoneutron multiplicities 
from the experimentally obtained average excitation energies has been made at 
160 < < 250 MeV in a model-dependent way [^ . 

For (7, n) channel the overall agreement between Livermore and Saclay data 
is good. Some inconsistency exists only in the giant resonance peak height (~ 3% 
for ^^'^Au and ~ 20% for ^'^^Pb) and on the right side of the peak, where (7, n) and 
(7,2n) process compete with each other, see Figs. ^ and ||. Large discrepancies 
are present in (7,2n) cross sections measured in different experiments, both in 
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shape and normalization, up to ~ 50% for ^''^Pb target. 

Several attempts of data evaluation and re-measurements have been made to 
reduce these discrepancies. Based on the observation that the total photoneutron 
yields, (7(7, n) + 2cr(7, 2n) + 3cr(7, 3n), obtained in Livermore and Saclay experi- 
ments agree well, an explanation for the discrepancy was put forward in Ref. [ pO[| . 
It was attributed to the different neutron multiplicity sorting procedures adopted 
in different laboratories. As was concluded in Ref. [^, the neutron multiplicity 
sorting procedure adopted at Saclay was not correct since some of the (7, 2n) 
events were interpreted as pairs of (7, n) events. 

In 1987 new measurements were made in Livermore ||31[, where it was found 
that the previously reported Livermore ^6] and Saclay ||23| results have to be 
rescaled. As was recommended in Ref. 
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m for i97Au 



the Saclay data of Ref. for both 
^^''Au and ^°^Pb nuclei have to be used with the correction factor of 0.93. We 
follow this prescription in using photonuclear cross section data in our RELDIS 
code. However, such a correction is not widely accepted and the authors of 
Refs. |TU|, |n|] use uncorrected Saclay data. 

Some new (7, 2n) cross section data were obtained in Ref. 
and ^''^Pb. For ^^^Au nucleus, (7, 2n) data are nearly the same as the Saclay 
(7, 2n) data of Ref. and the conclusion of Ref. concerning the neutron 
multiplicity sorting procedure seems to be not confirmed. Unfortunately, the 
recent data of Ref. [^] are available only up to a few MeV above the (7, 2n) 
threshold and the findings of Ref. can not be completely ruled out. 

One of the most recent measurements of (7, n) cross section for photoabsorp- 
tion on ^osp]^ ^^^g performed in Russia at the Saratov University [^. A fine 
structure of the low-energy wing of the Giant Dipole Resonance (GDR) has been 
investigated in detail. The photoneutron cross sections were obtained from the 
photoneutron yield curves by means of the statistical regularisation method. 

An evaluation of (7, n) cross section for ^osp]^ j^g^g i^qq^i made at the Moscow 
State University |Q by applying a statistical reduction method. Because of 
systematic uncertainties in calibration and normalisation, the general character- 
istics of the measured (7, n) cross section (the energy integrated cross sections, 
weighted-mean values) are different in different measurements ^ |26[. In the 
reduction method of Ref. the renormalisation corrections were introduced for 
both the energy and the cross section scales in order to obtain the best agreement 
between the general characteristics of (7, n) cross section measured in different 
experiments. 

The data obtained in Refs. [^, |33| for (7, n) reaction on ^"^Pb target are also 
plotted in Fig. ||. Good agreement with re-scaled Saclay data of Ref. |2^ is found 
up to (7, 2n) threshold. Unfortunately, the (7, 2n) reaction was beyond the scope 
of investigations in Refs. 
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It is evident from our consideration that the calculations of electromagnetic 
dissociation of ultrarelativistic ^^^Au and ^osp]^ nuclei can not be based exclu- 
sively on the photoneutron cross sections measured in experiments with real pho- 
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tons. Additional information on photoneutron cross sections for the whole energy 
domain of equivalent photons {Emm ^ < E^ax) can be obtained by using 
theoretical models of photoabsorption. This is particularly indispensable for dis- 
sociation channels with emission of many neutrons (> 3) and charged particles, 
p, (i, t, a,... 



4.2 Evaluation of photoneutron cross sections by using 
the GNASH code 

One of the two photonuclear reaction models used in the present work is the 
GNASH code [Q]. It is very precise in describing low energy neutron emission 
data [35 1, however, it can be used only up to the pion production threshold, 
at E^ < 140 MeV. Within this model the photoabsorption process is modelled 
through the excitation of the giant dipole resonance (GDR) at lower energies and 
the quasideuteron (QD) mechanism at higher energies. 

The photoabsorption cross section in the whole energy range from the thresh- 



old for neutron emission up to 140 MeV is thus written in the form [p5| : 



(Ta{E^) = cTcDniEy) + aQoiEy) 



(23) 



where (Tgdr is given by a Lorentzian curve with parameters taken from GDR 



systematics and corrected according to Ref. ||3T|. The latter term, (Tqdi is re- 
lated by a Levinger-type model to the experimental deuteron photodisintegration 
cross section, ad 
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NZ 

aQD{E,) = L^<Ja{E,)F{E,) 



(24) 



where A^, Z and A are, respectively, the neutron, proton and mass number of 
the corresponding nucleus. The Levinger parameter, L, is equal to 6.5, and 
F{E^) is a Pauli-blocking factor, which reduces the free deuteron cross section, 
ad{E^), by taking into account the Pauli blocking of the excited neutron and 



proton in the nuclear medium. In Ref. F was derived in the form of a 
multidimensional integral, approximated in the energy range 20-140 MeV by a 
polynomial expression: 



FiE, 



8.3714 X 10 



-2 



9.8343 X IQ-^E^ 



+4.1222 X 10~*E, 



4 772 



(25) 



-3.4762 X IQ-^'E^ + 9.3537 x lO'^E^ . 



and by an exponential one outside the considered energy range: 



F{E, 




73.3/E^), 
24.2/E^), 



E^ < 20 MeV 
E^ > 140 MeV 
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Thus, F{Ey) tends to zero if goes to zero, and to unity if goes to infinity 
and is continuous with Eq.(|5D at 20 and 140 MeV |35|. 



Finally, the experimental deuteron photodisintegration cross section is given 
by a simple parameterisation: 

cTdiE^) = Q1.2{E^ - 2.224)3/2/^3 ^ (26) 

where E^ is expressed in MeV, as in the previous formulae, and ad in mb. 

Due to the correlation between intranuclear nucleons in the absorption on a 
quasideuteron pair, the initial particle-hole configuration is assumed to be 2plh 
rather than 2p2h, see Ref. [^. In the GNASH code the initial interaction charac- 
terized by the total cross section of Eqs. (p^)-p5|) is followed by the preequilibrium 
emission of fast nucleons described by the exciton model Finally, when 



the nuclear system comes to equilibrium, sequential evaporation of particles is 
considered within the Hauser-Feshbach formalism . 

GNASH code results for (7,n), (7, 2n) and (7, 3n) cross sections are shown 
in Figs. and § for ^^''Au and ^osp]^ nuclei, respectively. Calculations describe 
(7, n) and (7, 3n) data very well. Taking into account existing disagreements 
between the results of different measurements of (7, 2n) cross sections, one can 



conclude that the GNASH results fall in between the Saclay |23] and Livermore 



data |]26[ for ^^^Pb, and very close to Livermore data for ^^^Au, that seems to be 
satisfactory for both cases. Therefore, one can use in Eq. (|^ the photonuclear 
cross sections, aA^E^), and branching ratios, f{E^,i), calculated by the GNASH 
code to estimate the mutual electromagnetic dissociation cross section cr^^(i | j). 
The influence of a constraint, E^ < Emax = 140 MeV, will be discussed in Sec. |[ 

4.3 Neutron emission simulated by cascade and evapora- 
tion codes 

Branching ratios for neutron emission in photonuclear reactions, f{E^,i), can 
be calculated by means of the extended cascade-evaporation-fission-multifrag- 



mentation model of photonuclear reactions [0 in the whole range of equivalent 
photon energies. Some details of the calculation method as well as numerous 
comparisons with experimental data used for the model verification were given in 
Refs. H, ^ . Here we describe only the general calculation scheme along with 
the modifications and advancements made in the model since the time when the 
works 1^ ^ ^ have been published. 

In the RELDIS model the values of the total photoabsorption cross section to 
be used in Eq. are taken from corresponding approximations of experimental 
data. In the GDR region the Lorentz curve fits were used for this purpose with 



parameters from Ref. |27| corrected according to the prescription of Ref. |31], as 
described in Sec. [4.1| . Above the GDR region, where the quasideuteron absorption 
comes into play, the total cross section is taken from the quasideuteron model of 
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Ref. |2l: 



N 7 

aQo{E,) = k^ar\E,). (27) 



Here is the meson exchange part of the cross section for the deuteron pho- 
todisintegration, '^d — *• np, and k is an empirical constant p4 |. 

Above the pion production threshold, at > 140 MeV, a universal behaviour 



(7A{Ey) oc y4 is observed (see Ref. |38] for the latest experimental data). This 
means that the total photoabsorption cross section per bound nucleon aA{E^)/A 
has almost the same energy dependence for light, medium- weight and heavy 
nuclei, C, Al, Cu, Sn, Pb, at least up to E^ 3 GeV. Therefore, having the data 
for one nucleus one can calculate the cross section for other nuclei. However, in 
this energy region the universal curve aA{E../)/A is very different from the values 
extrapolated from the cross sections on free nucleons, {Za^p + Na^n)/^, which 
are deduced from proton and deuteron data [Q. At E"^ > 3 GeV the universal 
behaviour breaks down, and the ratio aA{E^)/A for lead is 20-25% lower than 
for carbon |]39| , |40| due to the nuclear shadowing effect. In order to approximate 
the total photonuclear cross sections at E^ > 3 GeV we used recent results 
obtained within the framework of the Glauber-Gribov scattering theory and the 
Generalized Vector Dominance model ||3^, 0. Such calculations describe well the 
general trend of experimental data obtained for high energy photon absorption, 
although the data have very large uncertainties at E^ > 10 GeV. 

By comparing Tab. II of Ref. and Tab. |I| of the present paper one can find 
minor differences in the total ED cross sections due to using different parameteri- 
sations of aA{E^)/A at E^ > 3 GeV. Compared with the total ED cross sections, 
the single or double neutron emission cross sections are even less affected by the 
choice of the parameterisation. 

The RELDIS code performs the Monte Carlo simulation of the mutual disso- 
ciation process according to the following steps. First, a pair of the energies Ei 
and E2 of the photons exchanged between the colliding nuclei is selected accord- 
ing to the spectral function Afm{Ei, E2). Second, the photoabsorption process 
is generated in both nuclei leading to the formation of excited residual nuclei. 
Third, the de-excitation of both of the thermalized residual nuclei is simulated 
according to the statistical evaporation-fission-multifragmentation model, (the 
SMM model) [|. 

The evaporation of neutrons from excited residual compound-like nucleus is 
the main process responsible for the (7,?^-), (7,2n), (7, Sn) channels of photoab- 
sorption. The quality of description of such channels is very important for precise 
calculations of neutron emission in the mutual dissociation process. In the present 
paper the standard Weisskopf evaporation scheme is used with several mod- 
ifications taking into account the microscopic effects of nuclear structure in the 
nuclear mass and level density formulae. Such effects reveal themselves in the 
noticeable difference, up to ~ 10 — 15 MeV for heavy closed-shell nuclei, between 
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the values of the nuclear mass measured in experiments and those predicted by 
the macroscopic liquid-drop model. 

Moreover, this difference in mass, the so-called shell correction, and the level 
density parameter are strongly correlated. For closed-shell nuclei the actual val- 
ues of the level density parameter are essentially lower than the average values 
expressed as A/8 — A/ 10 MeV~^, and these values depend noticeably on the ex- 
citation energy. Proper accounting for these effects, as well as pairing effects, is 
important at low excitations, E* ~ 10 MeV, i.e. in the region where In and 2n 
photoemission processes occur. 

The above-mentioned shell effects are most pronounced at low excitation en- 
ergies, but almost disappear at E* > 50 MeV, see Ref. ||2^ for details. Several 
phenomenological systematics of the level density parameter were proposed to 
account for such behaviour, Refs. 0, Our calculations are based on re- 

sults of Ref. where all existing data on the level densities, decay widths and 
lifetimes of excited nuclei have been analyzed in the framework of the statistical 
model. 

However, the creation and subsequent decay of an excited compound nucleus 
formed after the photoabsorption in GDR region is not the only process responsi- 
ble for the neutron emission. Indeed, a giant resonance is a coherent superposition 
of (one-particle-one-hole) Iplh excitations. A particle or a hole can interact with 
another nucleon and create a 2p2h state. Further spreading to 3p3h states etc. 
finally leads to a statistically equilibrated system, the compound nucleus. In- 
stead of such evolution to equilibrium, a collective Iplh state can decay directly 
by the emission of one nucleon leading to a low-lying hole state in the residual 
nucleus, see among others Ref. [Q. After such direct In emission, the emission 
of a second neutron is generally impossible, even though the initial photon energy 
exceeds the 2n emission threshold. In such a way (7, 2n) channel is suppressed 
in comparison with the pure statistical decay. 

Although the GDR state in ^osp]^ nucleus decays mainly statistically, the 
existence of direct neutron emission has been clearly demonstrated in Ref. [^. 
The process where a fast nucleon is emitted and the final state ^°^Pb nucleus is 
left with low excitation energy < 3 MeV was identified in this experiment . 
An evidence of direct neutron emission in photoabsorption on Au and Pb nuclei 
was given in Ref. [^] based on the analysis of competition between In and 
2n emission channels. This was an independent confirmation of the findings of 
earlier works and devoted to the measurements of the neutron spectra in 
photoabsorption on the same nuclei. An excess of fast neutrons (kinetic energy 
> 4 MeV) with respect to the predictions of the statistical evaporation model 
has been demonstrated and attributed to the direct emission. 

In Ref. [44] a non-statistical contribution in excited ^osp]^ nucleus with 10 < 
E* < 30 MeV was successfully extracted. Out of this region the non-statistical 
contribution was found to be negligible. In our calculations we used the total 
fractions of the non-statistical neutron emission from Au and Pb nuclei as P^**" = 
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0.31 and 0.26, respectively, evaluated from experimental data in Ref. |]23|. Such 
values are in line with modern theoretical expectations |^ that the ratio of 
intensities of the direct and statistical neutron emission from photoexcited GDR 
in 208p]^ nucleus is about ~ 0.1. In the RELDIS code the emission angles B of 
non-statistical fast photoneutrons are generated according to the approximation 
VF(G) = A + B siia^ 6 which is found in Ref. We assumed that the direct 
In emission takes place at 7 < -E* < 22 MeV. 

Since the adopted P^^^ values have some uncertainties, we have investigated 
the sensitivity of results to these values. A part of calculations was made with 
P^*'' = 0, i.e. without accounting for direct emission, see Figs. |^, || and Tab. ^ 
As shown in Fig. ^ the (7, 2n) cross sections on gold calculated by the RELDIS 
code with P^^*" = 0.31 are very close to Saclay measurements while Liver- 
more results 1 26] are better described with P^**" = 0. Therefore the difference in 
calculation results obtained with P^*^ = and P^**" = 0.31 reflects the level of 
experimental uncertainties. 



4.4 Choice of cutoff parameter and nuclear density distri- 
butions 

Since the nucleon-nucleon interaction has an isovector component, the interfer- 
ence of nuclear and electromagnetic amplitudes cannot be excluded. Such in- 
terference term was considered in Ref. |^ and found to be small. Even for 
^^^Au nucleus colliding with heavy targets the interference correction to the sin- 
gle neutron removal cross section was found to be less than 0.5-0.6 % of the 
corresponding nuclear or electromagnetic contributions. Following this result, 
nuclear and electromagnetic parts of the dissociation cross section may be safely 
treated separately. In other words, one can add probabilities instead of coherent 
summation of amplitudes. 

Let us consider the way how the probabilities of the nuclear and electromag- 
netic contributions should be added to obtain the total dissociation probability. 
At grazing impact parameters relativistic nuclei are partly transparent to each 
other. Hard A^A^ collisions may be absent at all in the case of a peripheral event 
with a weak overlap of diffuse nuclear surfaces, while the electromagnetic inter- 
action may take place in this event leading to the electromagnetic dissociation. 
Generally, at a grazing collision either the nuclear, or electromagnetic interac- 
tion, or even both of them may occur. As an example of the latter case, a single 
neutron-neutron collision in the participant zone may lead to the neutron re- 
moval, while a photon may be emitted and absorbed by charged spectators in 
the same event. 

Therefore, in a detailed theoretical model a smooth transition from purely 
nuclear collisions at 6 <C Pi + P2 to electromagnetic collisions at 6 Pi + P2 
should take place. Such kind of transition was considered in a "soft-sphere" model 
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of Ref . . A similar approach was adopted in Ref . |]TT| , where the cross section 
for at least one type of dissociation, either nuclear, or electromagnetic, or both, 
was written as: 

a = 27r + V^^{b) - , (28) 

where ^""^(6) and V^^{b) are, respectively, the probabilities of the nuclear and 
electromagnetic dissociation at given impact parameter b. Putting explicitly the 
integration limits for each term one obtains: 



a = 27ry^ bdbV'"'\b) + 2'K J^^^bdbV^^{b) -2tt J^^^ bdbV'''"'{b)V^^{b). (29) 

Here the impact parameter cutoff values, 6""*^ and bf^, were used for the nuclear 
and electromagnetic interactions, respectively. However, due to several reasons a 
more simple expression is widely used P, ffl: 



a = a 



+ a""" = 27r / 6d6P""^(6) + 2n / M6P^^(6), (30) 

Jo Jbc 



where a single cutoff parameter, be, was chosen as: b^^ < be < ft""^- "^^^ ^^^^ 
reason is that with the latter condition one can effectively reduce the first and sec- 
ond terms of Eq. ( ^91) without subtracting the third nuclear-plus-electromagnetic 
term. Numerical results based on Eqs. (|29| ) and (|30| ) become similar as it was 



found for "sharp-cutoff" and "soft-spheres" models of Ref. |^ . Second, for heavy 
nuclei the difference between reasonable values of b^^, fo""'^ and b^. turns out to be 
less than 1 fm. As a result, the last nuclear-plus-electromagnetic term of Eq. (p9|) 
turns out to be small. Third, with Eq. ( PP] ) one can study the nuclear and elec- 
tromagnetic contributions separately. Therefore, independent parameterisations 
may be found in experiments for the nuclear and electromagnetic parts, o"""*^ 
and a^^. This is useful for studying nuclear and electromagnetic dissociation 
at ultrarelativistic colliders, where the products of nuclear and electromagnetic 
interactions populate very different rapidity regions, namely the central rapidity 
region and close to the beam rapidity, respectively. 

In the present paper the traditional form given by Eq. ( |5(]| ) is adopted with a 
common impact parameter cutoff, b^, for nuclear and electromagnetic contribu- 
tions. At relativistic energies, according to the widely used BCV parameterisation 



of Ref. [Q, be is estimated as: 

be = RBcviAY' + - Xscv{A\^"' + A-^^l^)). (31) 

The values Rbcv = 1-34 fm and Xbcv = 0.75 were found from a fit to Glauber- 
type calculations of the total nuclear reaction cross sections . 



The evidences in favour of the be choice according to the BCV parameterisa- 
tion were given in Refs. |]5D[,[^, which we mention among others. As argued in 



19 



Ref . , using the BCV parameterisation one can perfectly describe experimen- 
tal data on fragment angular distributions which are very sensitive to be- 

For calculations within the abrasion model we used the following values for the 
total nucleon-nucleon cross section given in Ref. ||52|: aNN = 40, 50 and 90 mb at 
SPS, RHIC and LHC energies respectively. Some problems are connected with 
the choice of the nuclear density parameters, -Rq = x A^^^, and d = 0.54 fm. 
Sec. 1^. Only the nuclear charge distributions are measured in electron scattering 
experiments, while the neutron densities are available only from calculations. 
We used = 1.14 fm as an average value between the proton and neutron 



distributions similar to one used in Ref. 11 . The total nuclear reaction cross 



sections calculated with these r„ and d values within the abrasion model are in 



good agreement with the approximation of experimental data found in Ref. [ p2 

Numerical results showing the sensitivity of the nuclear and mutual elec- 
tromagnetic dissociation cross sections to the variations of the above-discussed 
parameters are given in Sec. |[ 



5 Comparison with CERN SPS data 

The calculated charge changing cross sections of the single dissociation of 158A 
GeV 208p]^ jQj^g g^j^g given in Fig. ^. The LO mechanism with single photon 
exchange. Fig. |l|, and the NLO mechanism with double photon exchange. Fig. ^ 
were taken into account. 

Each proton removal process can be accompanied by neutron loss as well. 
The calculated cross sections given for Z = 82 correspond to the interaction 
where only neutrons are emitted. However, experimental data are absent for 
this channel. As one can see, the electromagnetic contribution dominates for the 
processes with removal of one, two and three protons where good agreement with 



the experimental data of Ref. |Q is found. 

Another important check of the model becomes possible with recent experi- 
mental data for Au fragmentation by ultrarelativistic Pb ions [Q. In this case 



the neutron emission is investigated directly. The experimental exclusive cross 
sections for emission of one or two neutrons are compared with theory in Tab. ^ 
Also in this case the contributions from the LO and NLO processes. Figs. |l| and 
1^, were taken into account. 

The calculations were made with and without accounting for direct neutron 
emission process, i.e. with P^*'' = 0.31 and P^'^^ = 0, respectively. The results 
with P^**" = 0.31 are in better agreement with experiment and this is especially 
true for 2n emission channel. Therefore, we use this value in further calculations. 

The predictions of abrasion model are also in good agreement with data [54 . 
Therefore, our choice of the critical impact parameter, be, is justified by such a 
comparison. As was stressed in Sec. ^ the interaction of knocked out nucleons 
with spectators and spectator de-excitation process themselves were neglected in 
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this version of the abrasion model aimed at considering only few nucleon removal 
processes. Such good agreement indicates that a simple abrasion model proposed 
for describing the data at ~ 1-10 AGeV can be used successfully at much higher 
energies as well. 

6 Mutual dissociation of ^^'^Au and ^'^^Pb ions at 
RHIC and LHC 

On the basis of the successful verifications at lower energies described in Sec. |], 
the model can now be extrapolated to RHIC and LHC energies. In a collider, 
the mutual heavy-ion dissociation process takes place at the crossing point of two 
beams. Downstream from this point the dissociation products can be separated 
by the magnetic field according to their Zj A ratio. Protons and nuclear fragments 
move close to the beam trajectories, while free neutrons leave the beam pipe after 
a dipole magnet. 

At RHIC, Zero Degree (very forward) Calorimeters (ZDC) for each beam are 
located after the magnets and they are well designed for the neutron registra- 
tion [P, |l^, [ill]. Therefore, in the following we consider mainly semi-inclusive 



mutual neutron emission cross sections, (Tm(^ | j), where i and j denote corre- 
sponding channels, InX, 2nX, 3nX,... Besides the emission of a given number 
of neutrons, such dissociation channels contain any number of other particles de- 
noted as X or y. Protons can be found most often amongst the particles emitted 
along with neutrons. The proton emission rates predicted by the RELDIS code 
were found to be in agreement with the data on Pb dissociation ||5^, see Sec. |[ 
Therefore we believe that our model is accurate in estimating Om{}-nX \ InF), 
ami^nX I 2nY) and am{'2nX \ 2nY) cross sections. 

If the dissociation channel of one of the collision partners in mutual disso- 
ciation is not exactly known, one can define inclusive mutual dissociation cross 
sections, am{lnX \ V), ami^nX \ V), (Tm(3nX | V) etc. In such notations V 
denotes an arbitrary dissociation mode of one of the nuclei. 

The cross sections for some specific channels of mutual dissociation at RHIC 
energies are given in Tab. |^. Compared to the total cross sections, the relative 
contributions of NLO processes are very different for the cross sections of In, 2n 
and 3n emission. As one can note, a^{lnX \ InY) has a small NLO correction, 
~ 7%, while (T^^(3nX | D) becomes almost twice as large as its LO value if NLO 
correction is included. This is due to the fact that the NLO processes shown in 
Figs. ^ and || include nuclear excitation due to double photon absorption and, 
particularly, the double GDR excitation process. Since the average GDR energy 
for Au and Pb nuclei is about 13-14 MeV, the double GDR excitation introduces, 
on average, 26-28 MeV excitation which is above the 3?7. emission threshold. 

Therefore, as a rule. In and 2n emission cross sections are less affected by 
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including of NLO corrections than 3n cross sections. The double GDR excitation 
process is a well-known phenomenon studied in low and intermediate energy 
heavy ion collisions, see, for example, Ref. [^]. However, nothing is known about 
the multiple photon excitations of heavy ions at ultrarelativistic energies, where 
the key question is to what extent the distribution of probabilities for multiple 
excitations follows the Poisson distribution (the harmonic picture of excitations). 

In order to estimate the rates of multiple excitations, including those well 
above the giant resonance region, the ratios a^^{2nX \ D)/a^^{lnX \ InY) 
and (T^^(3nX | D)/a^^{lnX \ InY) should be measured at RHIC. Such mea- 
surements do not require preliminary determination of the collider luminosity. 
As it follows from Tab. ^ the above-mentioned ratios are remarkably reduced 
if the multiple excitations, and hence NLO processes, are suppressed in mutual 
electromagnetic dissociation at RHIC. 

There remains still some freedom in choosing several parameters of our model. 
To check the sensitivity of our predictions to their variations we performed the 
calculations for a reasonable span of input parameters. 

Tab. ^demonstrates the sensitivity of the mutual electromagnetic dissociation 
cross sections to the photonuclear cross sections used as the input. In order to 
demonstrate such sensitivity we used two different models to calculate such cross 
sections, namely the GNASH code |^ and the photonuclear reaction model 
implemented in the RELDIS code |3^ itself. Additionally, in the latter model 
we used two different values for the probability of the direct neutron emission in 
In channel, -P^**". 

Besides the variations of the cross sections for emission of one or two neu- 
trons, the variations of a cumulative value, the Low Multiplicity Neutron (LMN) 
emission cross section defined as 



am{LMN) = a^ilnX \ InY) + a^{lnX \ 2nY) + 
a.^{2nX I InY) + am{2nX \ 2nY) 

were evaluated for several sub-regions of equivalent photon energies, < 24 MeV, 
< 140 MeV and for the full range. In the latter case the NLO processes were 
taken into account along with the LO mutual dissociation process. Therefore, 
the importance of NLO corrections can be also estimated from Tab. ^. 

By examining Tab. ^, one can draw several conclusions. First, the semi- 
inclusive cross section a^{lnX \ InY) = 437 mb calculated for the pho- 
toabsorption in the giant resonance region is very close to the exclusive value 
a^{ln I In) = 445 mb obtained in Ref. fill with the same condition: E^ < 24 
MeV. Second, the calculations based on the photonuclear cross sections predicted 
by the GHASH code are very close to the RELDIS results for the photon en- 
ergy region E^ < 140 MeV. The difference between the RELDIS results for 
E^ < 24 MeV and E^ < 140 MeV is explained by the contribution of the 
quasideuteron photoabsorption mechanism to the In and 2n emission. Third, 
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the calculations which take into account the quasideuteron photoabsorption and 
photoreactions above the pion production threshold give about 25% enhancement 
in a^^{lnX \ InY) if the whole energy region of equivalent photons is considered 
and the NLO corrections are properly taken into account. At the same time, the 
cross sections a^^{lnX \ 2nY) and a^{2nX \ 2nY) increase up to two and four 
times, respectively, compared with the values calculated at the GDR region. 

The cross sections for these dissociation channels are large and such channels 
can be easily measured in experiments. Although the LO process of photoab- 
sorption in the GDR region gives an important contribution, the whole range of 
the equivalent photon energies should be considered. Also the contributions from 
the NLO processes should be taken into account to obtain the realistic values of 
the dissociation cross sections. 

One more conclusion follows from the results presented in Tab. ^ Calculations 
with P^^^ = and P^^^ = 0.31 give 10-40% difference in specific dissociation 
cross sections, but the values of the LMN cross section, a^^{LMN), practically 
coincide. This cross section is very high, a^^{LMN) ~ 1100 mb, and thus can 
be used for luminosity monitoring. As was shown above in Sec. [4.1| , an inevitable 
systematic error of ~ 5% should be assigned to this value due to uncertainties in 
the photoneutron cross sections measured in experiments. 

Tab. ^ shows the sensitivity of the mutual dissociation cross section in grazing 
nuclear collisions to the variations of the Tq parameter in the nuclear density dis- 
tribution, Eq. ([T9|). The parameters of neutron density distributions are not well 
determined and this table demonstrates possible ambiguities in nuclear dissocia- 
tion cross sections caused by this fact. The cross section variations are smaller 
in Tab. |^ compared to Tab. ^ about 3-8%. A small decrease in the To parameter 
leads to a decrease in correlated In — In emission, but, on the contrary, it leads 
to an increase in In — 2n and 2n — 2n emission. However, the LNM cross section, 
a'^'^{LMN), turns out to be more stable, within ~ 2% variations, compared to 
the cross sections for specific channels. 

The sensitivity of the dissociation cross section in grazing nuclear collisions 
to the variations of the total nucleon-nucleon cross section is investigated in 
Tab. ^. The variations of the (Ti^n in the range of 40-60 mb have only a small 
effect, within 4%, on the cross sections of the specific neutron emission channels. 
Since such variations have different signs, the influence on the cumulative value, 
a'^^^LMN), is less noticeable, below 2%. 

The cross sections given in Tabs. ^ and § for grazing nuclear collisions were 
found to be lower compared with the electromagnetic dissociation cross sections 
of Tab. ^. Moreover, it should be stressed, that only a part of the neutrons in 
grazing nuclear collisions is emitted at the forward or backward angles covered 
by the ZDCs. Therefore, an exact relation between nuclear and electromagnetic 
dissociation channels in each heavy-ion experiment should be only obtained by 
taking into account detection acceptances and trigger conditions of the corre- 
sponding experimental setup. 
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Tab. 1^ demonstrates the sensitivity of the mutual dissociation cross section 
in electromagnetic and nuclear interactions to the impact parameter cutoff, be- 
This is an important input parameter which has a noticeable influence on the 
final result. By changing this parameter by 5%, within the range of 14.5-16 
fm, one obtains the variations of the electromagnetic dissociation cross sections 
within 3-8%. Such variations of be shift the point, which delimits the regions of 
nuclear and electromagnetic interactions, below or above the domain where the 
overlap of diffuse nuclear boundaries takes place. In other words, assuming first 
be ~ be^ and then be ~ fe""'^ and considering the difference in final results, one 
can prove the possibility to use Eq. (0) instead of Eq. p9| ) . 

For example, if the cutoff value, be, becomes lower, all the ED cross sec- 
tions. In — In, In — 2n, 2n — 2n and a^^{LMN) become lower. The variations 
of the nuclear cross sections are more noticeable and have the opposite trend: 
such cross sections become higher by 5-15%. Finally, as one can see in Tab. |^ 
am{LMN) variations for both types of interactions are weaker than the varia- 
tions of the corresponding individual cross sections, within 3-8%, while the sum 
ag^{LMN) + a'^^^LMN) is altered by 1-4% only. 

Concluding the investigation of the sensitivity of the final results to the model 
parameters, one can note that for AuAu collisions at RHIC a^^{LMN) = 1100 
mb and especially (j^'^{LMN) = 737 mb values are more stable with respect to 
such variations in comparison with the individual cross sections, am{^nX \ InY), 
ami^nX I 2nY) and others. 

The same tendency was found for PbPb collisions at LHC energies where such 
cross sections were found to be a^^{LMN) = 1378 mb and a'!^^{LMN) = 755 
mb. For both RHIC and LHC cases the overall uncertainty of the am{LMN) 
calculation method may be estimated at the level of 5-7%. 

The condition for heavy ion dissociation to be mutual leads to some specific 
features for nuclear and electromagnetic interactions. The former interaction 
causes mainly symmetric or quasi-symmetric dissociation. The latter makes very 
probable asymmetric dissociation like {InX \ 5nY) or even [InX \ lOnY). This 



feature is shown in Figs. 10 and O where the cross sections with one, two or 



three neutrons in one arm of the ZDC setup are presented. As one can see, 
{InX I lOnY) dissociation is almost absent in nuclear collisions. This is not true 
for electromagnetic dissociation where the number of {InX \ lOny) events is 
approximately 1-5% of the main dissociation channel, {InX \ InY). 

In the RELDIS model nuclei undergo dissociation independently in the elec- 
tromagnetic fields of each other. Therefore there is no correlation between the 
numbers of neutrons, ni and n2, emitted by each of the nuclei, and asymmetric 
dissociations are possible along with symmetric ones. The extreme case of asym- 
metric dissociation is, of course, the single dissociation process. The nuclear 
dissociation considered in the framework of the abrasion model has different 
characteristic features. Namely, the numbers of emitted neutrons and protons 
are correlated due to the condition zi + ni = 22+^^2, see Sec. |3|, and nuclear 
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dissociation is always mutual. The latter condition was also adopted in Ref. [p 



The main results of our study are presented in Figs. [T^ and |T3|. They show the 
electromagnetic and nuclear dissociation cross sections. Since In and 2n emission 
in electromagnetic collisions is enhanced due to the GDR and QD absorption 
mechanisms, the corresponding strips are prominent in the plots. Simultaneous 
GDR excitation in both of the nuclei is a dominant process leading to the mutual 
dissociation, but it is responsible only for a part of the total mutual dissociation 
cross section, ~ 18% at RHIC, for example. The rest is provided mainly by 
asymmetric processes, when one of the nuclei is excited in a GDR state, while 
another nucleus absorbs a photon with the energy above the GDR region, which 
leads to emission of many neutrons. 

As seen in Figs. |12] and |13|, the probabilities of the simultaneous emission of 
three and more neutrons are small and such processes with participation of high 
energy photons are distributed over the large area in the plots. 

For the sake of completeness the cross sections of the mutual dissociation with- 
out neutrons emitted from one or both of the collision partners are also included. 
The rates of such processes, when mainly the proton emission takes place, are 
small. This is another difference between the electromagnetic and nuclear disso- 
ciation. The nuclear interaction events, when only a proton is removed from one 
or both of the nuclei, are very probable, see Figs. and |13[ 



7 Conclusions 

Since its experimental discovery, the process of electromagnetic dissociation of 
heavy ions has been studied only in fixed target experiments. In the year 2000 the 
Relativistic Heavy Ion Collider (RHIC) became operational at BNL and among 
other experiments the electromagnetic dissociation of ultrarelativistic heavy ions 
can now be investigated in collider kinematics by means of the Zero Degree 
Calorimeters (ZDC) [^. This makes possible to study the mutual heavy ion 
dissociation process for the first time. 

In the present paper the equivalent photon method and the abrasion model 
for grazing nuclear collisions were extended to the case of mutual dissociation of 
collision partners. Our numerical results for the total heavy ion dissociation cross 
sections are very close to the results of related studies, Ref. 0, [llj. Restricting 
ourselves to the domain of equivalent photon energies below 24 MeV, where only 
the giant resonance excitations are possible in electromagnetic dissociation, we 



obtained the cross section of correlated In-ln emission close to one of Ref. [p 
However, as we have found, at collider energies the neutron emission process in 
mutual electromagnetic dissociation is not entirely exhausted by the simultaneous 
excitation and decay of the giant resonances in both of the colliding nuclei. Apart 
from the mutual GDR excitation, asymmetric processes with the GDR excitation 
in one of the nuclei accompanied by a photonuclear reaction in the other nucleus 
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are very probable. A wide set of photonuclear reactions should be taken into 
account to obtain a realistic estimation of the mutual dissociation rate. 

As we found, the cross sections of emission of two and especially three neutrons 
in mutual dissociation are very sensitive to the presence of the next-to-leading- 
order processes where a collision partner (or even both of the nuclei) absorbs 
a pair of photons in a collision event. The contribution of multiple photon ab- 
sorption processes to the total single electromagnetic dissociation cross section 
at ultrarelativistic energies is 1-2% only. This is very different from the case of 
mutual electromagnetic dissociation where 30-40% of events are due to multiple 
photon absorption. Indeed, the main part of mutual dissociation events takes 
place in close collisions with small b > where the probability to absorb a vir- 
tual photon is large, and hence, two or more photons can be absorbed by each of 
the collision partners. 

We have examined the reliability of our results by studying their sensitivity 
to the variation of input data and parameters. Trying to answer a key question 
on whether the mutual dissociation cross section can be calculated with high 
accuracy, we have critically reviewed our model assumptions and the results of 
previous theoretical and experimental studies of photonuclear reactions and heavy 
ion dissociation processes. 

The ambiguity in the calculations of In — In correlated emission cross section 
alone, a{lnX \ InY), is found to be up to 15%. This is mainly due to the 
difference in the values of the photoneutron cross sections measured in different 
experiments. However, the ambiguity is lower, ~ 5 — 7%, if the sum of one and 
two neutron emission channels, a{lnX \ lnY) + 2a{lnX \ 2nY) + a{2nX \ 2nY), 
is considered. Therefore, it is a kind of cumulative neutron emission rate which 
should be used as the luminosity measure at colliders. 

We have found several distinctive features of the mutual electromagnetic dis- 
sociation process which are helpful for its experimental identification. Beside the 
enhancement of In and 2n emission channels, the electromagnetic interaction 
leads to very asymmetric mutual dissociation channels where only one neutron 
is lost by one collision partner while many neutrons are lost by another part- 
ner. Such dissociation pattern is very unlikely in grazing nuclear collisions with 
participation of the strong nuclear forces. 

The correlated emission of one or two neutrons in both the forward and back- 
ward directions without any additional particles in the mid-rapidity region can be 
used as a clear sign of the electromagnetic dissociation of ultrarelativistic heavy 
ions. The identification of such mutual dissociation events and counting their 
rates in both arms of ZDC along with the calculation results of the present pa- 
per can provide a basis for an absolute luminosity calibration at RHIC. Similar 
methods can be used for the ALICE heavy ion experiment [57| planned at the 
future Large Hadron Colhder (LHC) to be build at CERN. 
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Table 1: Total cross sections (barn) of single and mutual dissociation calculated 
by the RELDIS code, abrasion model and by other authors for AuAu and PbPb 
collisions at RHIC and LHC 





Disso- 
ciation 

process 


First 
order 


Second 
order 


All 

contri- 
butions 


65+65 AGeV 
AuAu at RHIC 


Single 
electromagnetic 


82 


1.78 


83.8 


Mutual 
electromagnetic 




2.5 


3.6 


Nuclear 




- 


7.29 


100+100 AGeV 
AuAu at RHIC 


Single 
electromagnetic 


93.2 


1.86 


95.1 
88 i 
95 


Mutual 
electromagnetic 


0.39 • 10-3 I 
0.49 ■ 10-3 [ji 


1 


2.6 


3.8 
3.9 |] 


1 




Nuclear 






7.29 
7.09 


1] 




2.75+2.75 ATeV 
PbPb at LHC 


Single 
electromagnetic 


212 


3 


215 
214 i 
220 [0 


Mutual 
electromagnetic 


0.43 ■ 10-3 
0.54-10-3 


i 

12] 


3.9 


6.2 
7.15 


1] 




Nuclear 






7.88 



31 



Table 2: Cross sections (barn) for ^^^Au dissociation induced by 158 AGeV Pb 
beams. Theoretical results are obtained by the RELDIS code and abrasion model. 
Experimental data are taken from Ref. [Q. RELDIS results without direct In 
emission are given in parentheses. 



Dissociation 
channel 


af^(z)+af^(z) 


o-""^(i) 


All contributions 


Experi- 
ment 


RELDIS 

code 


Experi- 
ment 


Abrasion 
model 


Experi- 
ment 


Theory 


i = In 
i97Au ^196 Au + n 


26.4 ±4.0 


26.96 
(25.09) 


0.3 ±0.1 


0.43 


26.7 ±4.0 


27.39 
(25.52) 


i = 2n 
197 Au -^195 Au + 2n 


4.6 ±0.7 


4.57 
(6.39) 


0.13 ±0.4 


0.13 


4.7 ±0.7 


4.70 

(6.52) 
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Table 3: Mutual electromagnetic dissociation cross sections for AuAu collisions at 
RHIC. X and Y denote any particle, except neutron. V means any dissociation 
channel. Calculation results are given (a) for the leading order contribution only, 
(b) for the sum of leading order and next-to-leading-order contributions. 





Cross 
section 
(mb) 


(a) 

T r\ 


(b) 

T r\ 1 i\TT r\ 




ag^ilnX 1 


InY) 


612 


659 


65+65 AGeV 




\V) 


1244 


1502 


AuAu at RHIC 




\V) 


330 


446 




a^^{3nX 


\V) 


148 


274 




1 


InY) 


607 


652 


100+100 AGeV 


a^^ilnX 


\V) 


1257 


1518 


AuAu at RHIC 


ag^{2nX 


\V) 


341 


461 




ag^iZnX 


\V) 


155 


284 
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Table 4: Sensitivity of the partial mutual electromagnetic dissociation cross sec- 
tions to the variation of probability of the direct neutron emission in In channel, 
P^^^, to the input photonuclear cross sections and to the corrections for next-to- 
leading-order processes. Results obtained with GNASH and RELDIS codes are 
given for 100-1-100 AGeV AuAu collisions. Recommended values are presented 
in boldface. Prediction of Ref. for a^^{lnX \ InY) is given for comparison. 





< 24 MeV 
LO 


E.^ < 140 MeV 
LO 


Full range of 
LO+NLO 


Cross section 
(mb) 


RELDIS 

T)dir Q 

n 


GNASH 


RELDIS 

pdir Q 

n 


RELDIS 

pdir Q 

n 


RELDIS 

pdir = 0.31 


ag^ilnX 1 InY) 


437 

445 |11 


430 


467 


549 


652 


a^{\nX 1 2ny) + 
(J^{^nX 1 \nY) 


205 


221 


262 


439 


388 


o^{2nX 1 InY) 


21 


28 


38 


87 


60 


al^{LMN) 


663 


679 


767 


1075 


1100 
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Table 5: Sensitivity of the mutual dissociation cross section in nuclear interac- 
tions to the variation of Tq parameter of the nuclear density distribution. Results 
of the abrasion model are given for 100+100 AGeV AuAu collisions. Recom- 
mended values are given in boldface. 



Cross section 


Abrasion model 


(mb) 


To = 1.09 
Ro = 6.34 fm 


To = 1.12 
Ro = 6.52 fm 


ro = 1.14 
Ro — 6.63 fm 


To = 1.16 

Ro = 6.75 fm 


a'i^^lnX 1 \nY) 


361 


364 


371 


382 


aZ\^nX 1 2nF) + 
aZ\^nX 1 InY) 


241 


232 


224 


226 


G'^\2nX 1 2nY) 


148 


147 


142 


139 


al^%LMN) 


750 


743 


737 


747 
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Table 6: Sensitivity of the mutual dissociation cross section in nuclear interac- 
tions to the variation of the total nucleon-nucleon cross section a^N- The results 
of the abrasion model are given for 100+100 AGeV AuAu collisions. Recom- 
mended values are given in boldface. 



Cross section 
(mb) 


Abrasion model 


cnn — 40 mb 


cnn = 50 mb 


crjvAT = 60 mb 




InY) 


370 


371 


374 




2nY) + 
InY) 


233 


224 


222 




2nY) 


148 


142 


138 




751 


737 


734 
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Table 7: Sensitivity of the mutual dissociation cross section in electromagnetic 
and nuclear interactions to the critical impact parameter be- The results of the 
RELDIS code and abrasion model are given for 100+100 AGeV AuAu collisions. 
Recommended values are given in boldface. 



Cross section 
(mb) 


Rbcv — 1-27 
be = 14.45 fm 


Rbcv = 1-34 
be = 15.25 fm 


Rbcv = 1-41 
be = 16.05 fm 


ag^{lnX 1 InY) 


677 


652 


629 


ag^{lnX 1 2nY) + 
crg^{2nX \ InY) 


417 


388 


374 


ag^{2nX 1 2nY) 


62 


60 


57 


a^^(LMN) 


1156 


1100 


1060 


a'^^ilnX 1 \nY) 


379 


371 


390 


al'^^XnX 1 2nr) + 
(7'^''{2nX 1 InY) 


240 


224 


259 


(7^"^(2nX 1 2nY) 


141 


142 


151 


a'^^LMN) 


760 


737 


800 


a^^{LMN)+ 
a'^^^'^iLMN) 


1916 


1837 


1860 
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Figure 1: Electromagnetic excitation of one of the colliding nuclei: first order 
process. Open and closed circles denote elastic and inelastic vertices, respectively. 




Figure 2: Mutual electromagnetic excitation of relativistic nuclei: first order 
process. Closed circle denotes inelastic vertex. 




Figure 3: Mutual electromagnetic excitation of relativistic nuclei: second order 
process. Open and closed circles denote elastic and inelastic vertices, respectively. 
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Figure 4: Electromagnetic excitation of a single nucleus: second order process. 
Open and closed circles denote elastic and inelastic vertices, respectively. 




Figure 5: Mutual electromagnetic excitation of relativistic nuclei: next-to- 
leading-order contribution with single and double photon exchange processes. 
Open and closed circles denote elastic and inelastic vertices, respectively. 




Figure 6: Mutual electromagnetic excitation of relativistic nuclei: next-to- 
leading-order contribution with two double photon exchange processes. Open 
and closed circles denote elastic and inelastic vertices, respectively. 
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Figure 7: Photoneutron cross sections for gold. Open and closed circles are, 
respectively, Saclay and Livermore data re-scaled according to Ref. . 
GNASH code results are presented by solid line. RELDIS results are given by 
dashed and dotted lines for variants with and without inclusion of the direct In 
emission, respectively. 
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Figure 8: Photoneutron cross sections for lead. Open and closed circles are, re 
spectively, Saclay and Livermore |p6[ data re-scaled according to Ref. 
Saratov data 



Crosses 

tions are the same as in Fig. |^. 



3^, squares - Moscow evaluated data |Q. Other nota- 
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Figure 9: Charge changing cross sections of 158 AGeV ^osp^ ions on Pb target. 
The dashed- and dotted-hne histograms are the RELDIS and abrasion model 
results for electromagnetic and nuclear contributions, respectively. The solid line 
histogram presents the sum of both contributions. Points are experimental data 
from Ref. [El. 
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Figure 10: Mutual electromagnetic, cr^ (i | j) (top), and nuclear, a 



j) (bot- 



torn), dissociation cross sections for neutron emission in 100+100 AGeV AuAu 
collisions at RHIC. The cross section values for j = InX, 2nX, 3nX are given by 
the solid, dashed and dotted histograms, respectively. 
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Figure 12: Mutual dissociation cross sections (mb) for neutron emission = 
OnX, InX, ...19nX) due to electromagnetic {p^{i \ j), top panel), and nuclear 
(cr™'^(i I j), bottom panel), dissociation in 100+100 AGeV AuAu collisions at 
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Figure 13: Same as in Fig. O but for 2.75+2.75 ATeV PbPb collisions at LHC. 
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